Single-qubit gates and measurements in the surface acoustic wave quantum computer 
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In the surface acoustic wave quantum computer, the spin state of an electron trapped in a moving 
quantum dot comprises the physical qubit of the scheme. Via detailed analytic and numerical 
modeling of the qubit dynamics, we discuss the effect of excitations into higher-energy orbital states 
of the quantum dot that occur when the qubits pass through magnetic fields. We describe how 
single-qubit quantum operations, such as single-qubit rotations and single-qubit measurements, can 
be performed using only localized static magnetic fields. The models provide useful parameter 
regimes to be explored experimentally when the requirements on semiconductor gate fabrication 
and the nanomagnetics technology are met in the future. 
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I. INTRODUCTION 

Quantum computation promises enormous technolog- 
ical advances in the field of information processing^'^'^ii 
and the quest for its realization has attracted many 
strong contenders in the field of physics and engineer- 
ing. This paper is concerned with a scheme for quan- 
tum computation put forward by Barnes, Shilton and 
Robinsonpi which falls into the semiconductor quan- 
tum dot categoryi&LSi^iiSiiiii^ The proposal for quantum 
computation is based on the results of ongoing exper- 
iments that have demonstrated the capture and trans- 
port of single electrons in moving quantum dots^^*i^ii^ 
The dots are formed when a surface acoustic wave (SAW) 
travels along the surface of a piezoelectric semiconductor 
containing a two-dimensional electron gas (2DEG). See 
Fig. ^for a schematic diagram of the device. When the 
SAW is made to pass through a constriction in the form 
of a quasi-one-dimensional channel (QIDC), the induced 
piezoelectric potential drags electrons into and along the 
QIDC. In certain parameter regimes the device trans- 
ports one electron per potential minimum of the SAW.^** 
The spin on the trapped electron represents the physical 
qubit. Quantum computation involves performing qubit 
operations on the trapped electrons as they move with 
the speed of the SAW. 

Many schemes for quantum computation, such as con- 
ventional quantum dotS)Sii& doped silicon^ supercon- 
ducting boxeSfiSi and ion traps^ involve static qubits. 
The surface acoustic wave quantum computer on the 
other hand is of the "flying qubit" type, which include 
linear optics schemes)^ some ion trap schemes with ion 
shuttling^^ and schemes based on coherent electron trans- 
port in quantum wiresi ^i^^i^^ All these have in common 
that the carriers of quantum information physically move 
through space during the computation. Flying qubits 
have the advantage of being able to distribute infor- 
mation quickly over large distances across the quan- 
tum circuit when dccoherence times are short and to 
interface with quantum memory registers at fixed loca- 
tions. Another advantage of the SAW quantum com- 
putation scheme is its ensemble nature. It intrinsically 




FIG. 1: (Color online) Schematic diagram of an experimen- 
tal device for producing quantized acoustoelectric currents 
through a narrow QIDC constriction. 



performs time-ensemble computation, in much the same 
way NMR quantum computation performs molecular- 
ensemble computations i2^iSS*2& Time-ensemble computa- 
tion alleviates the demand for single-shot spin measure- 
ments and has the advantage of being robust against 
small random errors. 

This paper considers proposals^ for implementing 
quantum gates on single SAW qubits using only static 
magnetic fields generated by surface magnetic gates. De- 
tailed modeling of the gate operation has been accom- 
plished by means of both analytic solutions and numeri- 
cal simulations of the Pauli equation. We show how elec- 
trostatically confined moving electrons behave under the 
influence of various magnetic fields and discuss the im- 
plications for quantum computing with surface acoustic 
wave electrons. 

The physics of the SAW-guided qubit is explained in 
more detail in section ^1 In sections IIIII and IIVI we 
present results on single-qubit unitary gates and single- 
qubit readout gates based on the Stern-Gerlach effect. 
Section discusses some of the dccoherence processes 
involved in quantum-dot based schemes. Section IVII is a 
summary of the results with parameter regimes of inter- 
est for future experiments. 



II. SCHEME FOR QUANTUM COMPUTATION 

We begin by summarizing the quantum computation 
scheme proposed by Barnes, Shilton and Robinsoni^ Fig- 
ure n shows a schematic diagram of the experimental 
setup originaUy designed to demonstrate quantized cur- 
rents in semiconductorsiiii^ii^ A NiCr/Al interdigitated 
transducer is patterned on a GaAs/AlGaAs heterostruc- 
ture. A narrow depleted QIDC splits the 2DEG into two 
regions, the source and drain. When a high frequency 
AC signal is applied to the transducer, a SAW propagates 
through the 2DEG, producing a periodic piezoelectric po- 
tential across the 2DEG. The potential drags electrons in 
the source region through the narrow QIDC constriction 
into the drain. It has been shown experimentally that 
over a range of SAW power and gate voltages, the cur- 
rent passing through the QIDC is quantized in units of 
e/, where e is the electronic charge and / is the frequency 
of the SAWpi2J4,lft The smallest quantized current ob- 
served corresponds to the transport of a single electron 
in each SAW minimum. Typically, the SAW in GaAs 
moves at 2700 ms^^ at a frequency of around 2.7 GHz, 
with an applied power of 3 — 7 dBm4^ These parameters 
produce currents in the range of nanoamps. 

Given the ability to trap single electrons in the SAW 
minima, the scheme for quantum computation is as fol- 
lows. It is possible for an array of N QlDCs in paral- 
lel to capture TV qubits in every Mth minimum, with a 
single electron in each QIDC, producing a qubit regis- 
ter along the SAW wavefront. M can be chosen suffi- 
ciently large to ensure that the Coulomb interaction be- 
tween successive qubit registers do not interfere with each 
separate computation. The qubits move with the min- 
ima of the SAW, passing through a sequence of static 
one- and two-qubit gates before arriving at an array of 
spin readout devices. One-qubit gates may be operated 
by nanoscale electromagnetic fields. Where two-qubit 
gates are needed, neighboring QlDCs are allowed into 
a tunnel contact controlled by a potential on a surface 
gatei^ The use of the Coulomb coupling between neigh- 
boring QlDCs is a common tool in spintronics which can 
be used to generate entangled states in dual-rail qubit 
representations iSLSS Figure |21 illustrates the network of 
QlDCs and qubit gates envisaged for performing a par- 
ticular quantum computation. 

The SAW-trapped electron is well confined in all three 
spatial directions. In this paper we define Cartesian co- 
ordinates such that the SAW propagates along the x axis 
with the z axis normal to the 2DEG. The 2DEG is pro- 
duced by a band-energy mismatch at the GaAs/AlGaAs 
interface which gives rise to a confining potential in the z 
direction. The energy level spacings in the 2DEG well are 
on the order of 50 — 100 meV3* Further confinement in 
the y direction is provided by an extended QIDC, etched 
into the surface so as to avoid screening the SAW-induced 
potential with metallic surface gatesi^ Finally, the con- 
finement in the x direction is due to the SAW poten- 
tial minimum which is approximately sinusoidal. The 
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FIG. 2: (Color online) Schematic diagram of a quantum gate 
network in a SAW quantum computer: (a) Two-qubit tun- 
neling gates; (b) One-qubit magnetic gates (in various ori- 
entations); (c) Gate network for quantum computation with 
SAW electron spins. Gray lines running horizontally repre- 
sent QlDCs; blackened regions indicate the SAW minima 
where the qubits reside; arrows represent spin polarization; 
rings represent magnetic surface gates; white squares repre- 
sent readout gates. 



SAW amplitude is typically 40 meVpi which is suffi- 
ciently large to prevent qubits being lost via tunneling 
into neighboring SAW minima. 

There are two important aspects of the SAW quan- 
tum computation scheme that distinguishes it from other 
similar quantum dot schemes. First, the scheme pro- 
vides repetitions of the same quantum computation with 
each passing of a single wavefront of the SAW. There- 
fore, a statistical time-ensemble of identical computa- 
tions can be read out at the end of the QlDCs as a mea- 
surable current, alleviating the need for single-electron 
measurements. Two sources of noise in the measured 
current may be estimated as follows: The shot noise is 
largely determined by how well the current is quantized 
to J = e/sAW and precisions of < 0.1% can be exper- 
imentally achievedi^ Johnson noise arises from the re- 
sistance of the ohmic contacts and the 2DEG which are 
on the order of ^ 100 — 1000 51. At temperatures of 
1 K they produce a rms voltage noise spectral density on 
the order of 10~^" V/vHz at most. This cannot drive 
a current noise through the SAW device since its effec- 
tive intrinsic impedance of 10 Mf2 is comparatively very 
large. See experimental papersi^ti^ii^iSSiS for more de- 
tailed discussions. The second key aspect of the scheme 
is the static nature of the gate components of the quan- 
tum circuit. This alleviates the need for strong, targeted 



and carefully-timed electromagnetic pulses that can be 
difficult and expensive to implement. The requirement 
of such expensive control resources often limit the scala- 
bility of most quantum computing implementations. 

It would certainly be convenient, though not essential, 
to have a means of preparing a pure fiducial qubit state. 
In NMR quantum computing, operations are carried out 
on ensembles of replica qubits which remain close to a 
highly mixed state of thermal equilibrium. Nevertheless, 
a successful readout of the computation is obtained be- 
cause the sum over many identical computations provides 
a measurable signal. Similarly in the SAW quantum com- 
puting scheme, states close to the maximally-mixed state 
are still useful because of the time-ensemble nature of the 
scheme. However, in contrast to NMR schemes, it is in 
principle possible to read out single electrons in the SAW 
scheme. We will therefore begin with nonensemble quan- 
tum computation in mind and only later exploit the ad- 
vantages of ensemble computation to deal with noise. Of 
course, the more pure the qubit states remain, the faster 
the computation will converge to the result. For these 
reasons, we describe in this paper two simple methods 
for preparing pure fiducial qubit states. 

In the original proposal for the SAW quantum 
computer f^ it was noted that the application of an exter- 
nal magnetic field of about 1 T will influence which spin 
states of electrons are favored in the capture process. '^^ 
This polarized capture process can be summarized as fol- 
lows. The SAW is strongly screened in the bulk 2DEG 
until it is raised above the Fermi energy and the quan- 
tum dot begins to form. When the higher-energy minor- 
ity quantum dot forms for the higher-energy sub-band of 
polarized electrons, the probability of capturing minor- 
ity electrons is small (see Fig. |31), whilst the probability 
of capturing electrons from the lower-energy sub-band is 
large. Once a cloud of approximately polarized electrons 
is captured, the exchange interaction will generally en- 
tangle electrons in the same dot together, so that the 
subsequent loss of electrons into the Fermi sea will lead 
to a decoherent process that could relax the remaining 
electrons into the low-energy polarized state. A more de- 
tailed multiparticle analysis will be needed to determine 
the final state of the remaining single electron, but the 
combination of the above two processes is likely to lead 
to a high level of polarization. 

The above method is conceptually simple and would 
be easy to implement in the laboratory. However, the 
macroscopic magnetic field required in the capture re- 
gion will need to be shielded from the rest of the device 
where the quantum computation is to be carried out, and 
this may present a nontrivial problem. If we chose to 
drive spin flips using microwave pulses, then the macro- 
scopic field is actually required across the whole device. 
However, the problem with using microwaves is their rel- 
atively long wavelength which tends to affect every part 
the computation. Alternatively, we may use local, static 
magnetic fields to initialize, rotate and read out single 
qubits without a global magnetic field. In Section Hvl it 
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FIG. 3: (Color online) Energy diagram showing polarized 
electron capture by means of Zeeman band splitting in the 
presence of a uniform magnetic field: Arrows on electrons in- 
dicate spin polarizations; Ef is the Fermi energy. The SAW 
is strongly screened in the bulk 2DEG below the Fermi en- 
ergy. As the SAW enters the QIDC constriction, the confining 
potential begins to form. The probability of capturing spin- 
down is small at the point the minority spin-type dot forms 
(upper curve). In contrast, the probability of capturing spin- 
up electrons is high when the majority spin-type dot forms 
(lower curve). 



will be demonstrated that spin-polarized electrons can be 
prepared and measured using a gate driven by the Stern- 
Gerlach effect. First we turn to the implementation of 
single qubit rotations using local static magnetic gates. 



III. SINGLE-QUBIT UNITARY GATES 

Single-qubit unitary operations may be carried out, 
in principle, by allowing the trapped single electrons 
to pass through regions of uniform magnetic field. If 
there is no spin-orbit coupling effect, the spin state of 
the qubit will evolve according to the Zeeman term in 
the Hamiltonian: hgfiB B ■ cr ip, where B is the mag- 
netic field, cr is the vector form of the 3 Pauli oper- 
ators and ip = (a, /3) is a spinor. The Bloch vector 
n = (-2Im[a^*],2Re[a/3*], |ap - |/3p) precesses about 
the direction of the magnetic field with angular frequency 
g^B\B\/h. A local static magnet may be produced by 
a magnetic force microscope or by evaporative deposi- 
tion of a ferromagnetic material such as Cobalt, or a 
permalloy such as NiFe. Inevitably, different samples 
will produce different strengths of magnetic field. But 
there are methods to vary the strength and pattern of 
the field once the sample has been fabricated. It has been 
demonstrated that ferromagnetic properties of thin-film 
3d transition metals can be modified via ion irradiation^ 
Another method would be to use oxidation techniques 
with the atomic force microscopei^^*^Si2i Only two inde- 
pendent directions of the B-field are necessary to pro- 
duce an arbitrary single-qubit manipulation and we will 



choose these to be perpendicular to the direction of the 
SAW, one aligned with and the other perpendicular to 
the 2DEG. 

For idealized qubits with no spatial degree of freedom, 
the above model for single-qubit rotations is complete. 
However, the trapped electron is a charged particle with 
a spatial distribution within the dot. The fields couple 
to both the spatial and spin degrees of freedom, causing 
the electron to experience the Lorentz force as well as 
spin-precession. It is therefore clear that one cannot in- 
crease magnetic fields arbitrarily, since the Lorentz force 
will upset the confinement properties of the electron. Nor 
can the direction of the field be chosen arbitrarily without 
consequences for the robustness of the gate. To address 
these concerns, we analyze the behavior of the full wave- 
function: a two-component spinor field defined over space 
and time, il^a{x, t), under the action of gates operated by 
static magnets. 



A. Pauli Hamiltonian with uniform magnetic fields 

Assuming a uniform magnetic field in the lab frame, 
we solve the Pauli equation for the spin field from which 
the probability density field and the Bloch vector field 
can be obtained. 

The qubit is trapped in a net electrostatic potential 
with contributions from the QIDC split gates, the 2DEG 
confining potential and the SAW piezoelectric potential. 
The parameter regimes we consider allow us to neglect 
motion out of the 2DEG plane. ^^ The net confining po- 
tential in the xy plane is modeled by^l 



acting on the particle, the Schrodinger equation simpli- 
fies significantly and there would be no need to concern 
ourselves with the spatial behavior of the qubit. 

Before proceeding further, it is convenient to transform 
into the rest frame of the electron moving with the SAW 
speed V. We may use simple Galilean transformations 
x' = X — vt, y' = y, t' = t, since the SAW velocity is 
non-relativistic. The potentials in the electron rest frame 
become, in the harmonic oscillator approximation, 
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with k = 27r/A. A further convenience is to use a system 
of natural units such that h, m*, v, e are unity. The units 
of length, time and energy become h/m*v — 0.640 /im, 
h/m*v'^ — 0.237 ns and m*w^ = 2.78 ^eV, respectively. 
The natural unit of magnetic field is 1.61 mT. Parame- 
ters can now be assigned with dimensionless values with 
respect to the above units. In the electron rest frame. 
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where the primes will subsequently be dropped from the 
coordinates. We will further simplify the Hamiltonian by 
adopting the Coulomb gauge V • A = 0. 
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The QIDC split gate voltages are such that typically 
Vb -- 2800 meV. The width of the QIDC w is typically 
between 1 and 2 /xm. The amplitude of the SAW is A '--^ 
40 meV with wavelength A ^ 1 /im and frequency / ~ 
2.7 GHz. 

The appropriate nonrelativistic equation for the two- 
component spinor field is the Schrodinger equation with 
a Pauli Hamiltonian. For an electron moving in an ar- 
bitrary vector potential A and potential energy V, the 
Pauli Hamiltonian H is 
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where e > is the electronic charge, g ~ 0.44 is the 
Lande g-factor and m* ~ 0.067TOe is the effective mass 
in GaAs. By simulating in the plane of the 2DEG we sim- 
plify to a 2-1-1 dimensional model, with variables (x, y, t). 
If the particle were chargeless but with an anomalous 
magnetic moment, terms in H that explicitly involve e 
can be dropped. In such a case, there is no Lorentz force 



B. Qubit rotation: Uniform transverse magnetic 
field 



For a magnetic field in the y direction in the plane of 
the 2DEG and transverse to the QIDC, we model the 
magnetic field in the region of the single-qubit gate with 
a vector potential of the form 



A.{x) 



-j3x, 
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with Ax — Ay — 0. Clearly this does not vanish at in- 
finity, but we only consider interactions over regions of 
finite extent. This potential generates a uniform mag- 
netic field of strength /3 in the y direction. With this 
potential Eq. ((HJ becomes 

-^V27A + i/32(a:+t)2V'+igMB/3a,y7A-fF^ = z— tA, (8) 

in which the potential V is given by ^. An effective 
potential V ^^ A? j2 can be identified, which resembles a 
harmonic oscillator but is time-dependent. Anticipating 
an analytic solution by separation of variables we apply 
the ansatz 
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where 0„ are the harmonic oscillator eigenstates with 
energies E^ = ujy{n + ^), and \sy) are the eigenstates 
of ay with eigenvalues Sy = ±1. We have introduced 
the oscillator frequency tOy — ■\/Vo/u' and the Zeeman 
energy gap AE — g^L^p. The energy eigenstates in the y 
direction are exactly the harmonic oscillator modes. The 
wavefunction in the x direction is time-dependent and 
it is the solution of x(a;,t) to which we now turn. On 
substituting Q into (O, we obtain a PDE for x(a;7i): 



-\i>lx 



-(co + ci)x^ 



Ci xt 



X^idtx, (10) 



with Co = Ak^ and ci ^ 0^ . Terms which depend only 
on t have been dropped, as they merely contribute global, 
time-dependent phases that do not affect the dynamics. 
A Gaussian solution of H10() can be found with a further 
ansatz 



X{x, t) = exp [hit) x^ + /2(i) X + hit)] 



(11) 



This system of time-dependent functions can be deter- 
mined self-consistently, assuming a Gaussian groundstate 
of the SAW dot for the initial condition at i = 0. The 
resulting expressions for fi,hif3 are complicated, but 
by noting a few general properties of the solution we can 
understand all the important features of the dynamics. 
h{t) takes account of the normalization but is otherwise 
of no more interest. /i(t) and h{t) together describe a 
groundstate Gaussian wavefunction evolving in the time- 
dependent vector potential. The solution is particularly 
simple in that it remains Gaussian throughout, so we will 
only need to keep track of the position of the central peak 
/i and the width (standard deviation) S of the probability 
distribution |x(a;,t)p: 
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The typical energy scales for the SAW electron encoun- 
tered in experiment puts the system in the regime where 
s <C 1 and ci <C Co. Expanding /i and h to first order in 
s and ci/co, we obtain the asymptotic behavior for the 
position of the peak: 
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From this result we see that the magnetic field introduces 
a constant drift velocity of the peak in the —x direction. 
On exiting the interaction region, the peak will be off- 
center with respect to the SAW dot and the probability 
distribution will subsequently oscillate in the x direction, 
perhaps exciting higher-energy orbital states of the dot. 
If the charge distribution is pulled too far off-center, it 
is likely to escape the SAW quantum dot. We could ask 
how long the electron can remain in the field before it is 



dragged a distance A away from the center of the dot in 
the X direction: 
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When this is compared with the time required for the 
Bloch vector to rotate by some appreciable angle such as 
a TT rotation, we obtain the ratio 
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which is about 2000 for a 1 T field. This means that the 
Bloch vector is allowed to precess approximately 1000 
times before the drag effect destabilizes the qubit confine- 
ment. In the limit of weak fields, this ratio is no smaller 
than g/XB Ak ^ 800. Therefore, whilst it is possible that 
higher orbital states are excited as a consequence of the 
interaction, there is little danger of the electron escap- 
ing the SAW quantum dot during the operation of the 
gate. We now turn to the behavior of 6{t), the width of 
the Gaussian wavefunction, which in the limit ci <C cq 
oscillates according to the simple expression 
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with Cq — Ak"^ and ci = /5^ . The frequency of the oscilla- 
tion increases with the energy of the dot and is typically 
much faster than the Zeeman spin precession frequency, 
i.e., rate of precession of the qubit Bloch vector. These 
results are plotted for a specific case in Fig. 0] The top 
two plots show the evolution of n and 5. As the Bloch 
vector rotates about the field along the y-axis, the Gaus- 
sian translates in the x direction with a rapidly oscillating 
width. The bottom plot of Fig. 0] shows probability am- 
plitudes of excitation into higher-energy simple harmonic 
oscillator (SHO) modes in the x direction of the quantum 
dot. 

Our conclusions are as follows. In addition to rotat- 
ing the Bloch vector as required, the gate has the ef- 
fect of displacing the center of the Gaussian wavefunction 
which will increase the energy of the bound state by an 
amount depending on both the gate time and the field 
strength. Higher-energy orbital states are likely to be 
excited, which could lead to decoherence via spin-orbit 
couplings or dipole coupling to phonons and photons. 
According to Fig. 0] the electron is almost completely 
out of the groundstate after a 7r/2 rotation of the Bloch 
vector. In the extreme case after many spin precessions 
(typically thousands), the qubit will leave the SAW quan- 
tum dot. However, for quantum computations the gate 
time need only be long enough to conduct a single or- 
bit of the Bloch sphere, in which case the excitation into 
higher orbital states is negligible. These observations re- 
veal some important features of the gate operation that 
are not revealed by an idealized spin-only model of the 
qubit. 
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FIG. 4: Gaussian evolution of the probability distribution 
|x(a^,i)P in a constant magnetic field By. The top plot shows 
peak position /i and the middle plot shows the width S. The 
bottom plot shows the probability amplitudes \C„\^ of the nth 
SHO mode in the x direction: |Cip (long dash), |C2|^ (short 
dash). The plot also shows 1 — I Co I ^ (solid). Typical values for 
the SAW dot are taken: A = 40 meV; A = 1.0 ^m; magnetic 
field of 1 T. T — -n /g^sP ~ 0.08 ns is the time taken for a 
TT rotation of the Bloch vector about the y-axis. The initial 
Gaussian wavefunction was taken to be the SAW quantum 
dot groundstate (Co = !)• In the limit ci/co <C 1 (weak 
field), ^{t) moves approximately linearly and 5{t) undergoes 
bounded oscillations about 5(0) — s/-\/2. 



C. Qubit rotation: Uniform perpendicular 
magnetic field 

To move the qubit state to an arbitrary point on the 
Bloch sphere a second axis of rotation on the Bloch 
sphere is needed, and to this end we consider a gate im- 
plemented by a uniform magnetic field in the z direction. 
This will not be just a trivial extension of the preceding 
analysis, since the 3D rotation symmetry is broken by 
the motion of the SAW in the x direction. In practice, 
the y magnetic field may be easier to fabricate than the 
z magnetic field, since the latter passes perpendicularly 
through the 2DEG structure. It could feasibly be pro- 
duced by layering oppositely aligned thin-film magnets 
just beneath and just above the 2DEG, or alternatively 
by applying a global z magnetic field which is shielded in 
regions where it is not needed. 

A vector potential generating the uniform Bz magnetic 
field is 



^x(z/) = -/3y, 



(18) 



with Ay = Az = 0. It should be noted that although 
the field is uniform and static in the laboratory frame, 
the electron sees a moving uniform field. In the electron 
rest frame. Ax is time-independent and z-independent, 
allowing a Hamiltonian in x and y only. The chosen gauge 
facilitates numerical simulations which follow shortly. An 
effective scalar potential that is quadratic in y arises from 
the A^ term in the Hamiltonian (0. It is interesting to 
compare the strength of this confining potential with the 
QIDC potential, which is also approximately quadratic 
iny: 
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With magnetic fields of order 1 Tesla and typical QIDC 
energies, this ratio is of the order unity. This means that 
for the parameter values being considered the effective 
scalar potential arising form the A^ term is compara- 
ble to the QIDC confinement potential. A stronger field 
would start to significantly deform the shape of the dot. 
This will not be a problem as long as the evolution has 
occurred adiabatically during the deformation. However, 
as we shall discuss in a moment, the probability of exci- 
tation into higher orbital states due to the B^ field is not 
negligible. 

The other interesting term in the Hamiltonian is the 
asymmetric coupling —iPy dx between the x and y vari- 
ables. This is expected to introduce rotational behavior, 
as we would intuitively anticipate some form of Landau 
orbital motion due to the Lorcntz force. 

A numerical simulation implementing a Crank- 
NicholsonjS^ finite difference algorithm (alternating di- 
rection method) was used to simulate the operation of the 
gate. We started with an initial Gaussian groundstate of 
the dot with spin state | t^) and subjected it to 1 T of 
magnetic field for a duration of T^/2 = 7r/2g/iB/3, which 



is the gate time required for a 7r/2 rotation of the Bloch 
vector about the z axis. The evolution of the Bloch vec- 
tor is simple due to the uniformity of the magnetic field: 
nxii) — cos{g^iBPt), nyit) = sin((7/iB/3t), and tizlt) = 0. 
All other parameters were assigned those values given 
just after Eq. |^. The result of the simulation is shown 
in Fig. [SI which shows time-shots of the probability dis- 
tribution in the 2DEG {xy) plane. Losing its initial el- 
liptic contours, the distribution develops two lobes which 
rotate about its midpoint. The density at the center in- 
creases due to the spatial squeezing from the A^ term 
in the Hamiltonian. This clearly shows excitation into 
higher-energy orbital states. 

Using perturbation theory to second order in the 
strength of the field /?, with harmonic oscillator modes 
as the basis, we found the following facts: (i) To sec- 
ond order of perturbation theory, only the second excited 
state with SHO quantum numbers Ux — riy = 1 becomes 
populated; (ii) The probability ratio with respect to the 
groundstate is 



Cii 



Co 



= /3' 






(20) 



where tOx and ujy are frequencies arising from the har- 
monic oscillator approximation; (iii) The amplitude of 
the ratio approaches unity when the field approaches 1 T; 
(iv) The ratio of the spin precession frequency to the fre- 
quency of the |Cii| oscillation is about 1.6 x 10~^ for 
typical SAW parameters (see just after Eq. ^). 

The above model is useful in assessing the robustness 
of the qubit and its susceptibility to decoherence due to 
orbital motion. Population into this excited state is a 
problem for decoherence, since the oscillating charge in 
the dot couples via dipole interactions to phonons and 
other charges outside the dot. However, provided we 
have sufficient control of the gate time, we can use (|2Uf) 
to make the electron exit from the gate in its ground- 
state. Otherwise, a gate driven by radiofrequency pulses 
in the presence of a global magnetic field could provide 
an alternative means to implement spin rotations about 
the z axis. 
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FIG. 5: Contour plots of electron probability density during 
the operation of a qubit rotation gate. This is a numerical 
simulation in the xy (2DEG) plane of a SAW electron under- 
going a ■k/2 rotation gate about the z axis, under a uniform 
magnetic field B^ of strength 1 T. Snap-shots are shown at 
the following times (ps): (a) 0, (b) 2.89, (c) 11.6, (d) 20.3, 
(e) 28.9, (f) 37.6. The initial probability distribution is a 
Gaussian groundstate with standard widths a^ ~ 26.0 and 
ay — 20.6. The probability density rotates about the z axis 
under the infiuence of the Lorentz force acting in the xy plane. 
Moreover, an effective scalar potential contributed by the A^ 
term in the Hamiltonian is quadratic in y, which spatially 
squeezes the initial Gaussian distribution into the center of 
the dot. 



IV. SINGLE-QUBIT INITIALIZATION AND 

MEASUREMENTS 



In addition to single-qubit rotation gates, we require 
the ability to initialize and measure qubits at the be- 
ginning and end of the computation. Spin-polarized 
electrons can be obtained from injection through a 
ferromagnetic contact. '^^■'*° There is also a method to 
polarize spin using nondispersive phases (Aharonov- 
Bohm and Rashba) without the need for ferromag- 
netic contacts4i In the field of quantum computing, a 
well-known method for achieving readout of solid state 
spin qubits is to convert spin information into charge 
informatior^ and subsequently use single-electron tran- 
sistors or point contacts to detect charge displacements or 



Rabi oscillations 42il^ However, recent theoretical results 
by Stace and Barrett^ argue the absence of coherent os- 
cillations in a continuously measured current noise, con- 
trary to previous results and assumptions^^*^ and there- 
fore raise concerns about the measurability of charge os- 
cillations in similar scenarios. In any case, it is difficult to 
apply these methods to qubits in moving quantum dots. 
We therefore turn towards a quite different approach, 
one which enables the initialization and readout of SAW 
electron spin qubits solely with the aid of nanomagnets 
and ohmic contactsi^ The readout gate we consider is 
based on the Stern-Gerlach effect4L^2i^ In the 1920s 
Bohr and Pauli asserted that a Stern-Gerlach measure- 
ment on free electrons was impossiblej^Si^ using argu- 
ments which combined the concept of classical trajecto- 



ries and the uncertainty principle. This subsequently led 
physicists to analyze single-electron Stern-Gerlach mea- 
surements within increasingly more rigorous quantum 
settings, ultimately ending the debate by showing that 
the measurement can indeed be done, albeit with cer- 
tain caveatsi^SiSiSi Thus Stern-Gerlach measurements on 
free electrons have been extensively investigated, but lit- 
tle attention has been given to such measurements on 
confined electrons i^ An interesting semiclassical analy- 
sis of a Stern-Gerlach type experiment with conduction 
electrons has been reported;^ in which the authors jus- 
tifiably neglect the Lorentz force effects. In contrast, we 
analyze a single-electron Stern-Gerlach device, providing 
a full quantum mechanical treatment and including all 
Lorentz force effects. 

In the SAW electron system, single electrons are trans- 
ported in moving quantum dots. Can this system un- 
dergo a Stern-Gerlach type measurement, and would it 
be of any use in a scheme for quantum computation using 
qubits trapped in surface acoustic waves? 

The electron confinement to low dimensions allows us 
to guide the electron through the magnetic field in ways 
that enhance the spin measurement and suppress the 
deleterious effects of delocalization and the Lorentz force. 
Surface magnets can be arranged in such a way as to pro- 
duce a local magnetic field inhomogeneity. For example, 
two north poles placed on either side of the QIDC will 
produce a region of intense magnetic field gradient in- 
between (see Fig. Ejl. Via the Zeeman interaction term 
(X a ■ B in the Pauli Hamiltonian, the field inhomogene- 
ity has the effect of correlating the spatial location of a 
wavepacket to its spin statei^ In most situations it is 
necessary to continue confining the qubit during the op- 
eration of the gate, because the spreading time'^^ for a 
free Gaussian wavefunction is comparatively short - on 
the order of 1 ns. But the gate must cause a wavepacket 
splitting in order for the spin states to be resolved, hence 
the QIDC must relax to allow for motion in the y di- 
rection. This is achieved by patterning the QIDC in a 
funnel shape, with an opening angle 9 (see Fig.E)), such 
that in the electron rest frame the potential looks like 



^QlDC = Vo 2/V2 [W + tf,Il{e)u{t)tY 



(21) 



where u{t) is the step function: u{t < Q) — and 
u{t > 0) — 1. A negatively biased surface gate placed on 
the X axis can be used to guide the electron into the chan- 
nels. If necessary, both the position of the electrode and 
the opening angle 9 could be optimized for a particular 
sample device under low-temperature and high-vacuum 
conditions by erasable electrostatic lithographyiS We 
will study the quality of the readout obtained from the 
gate for two different magnetic field configurations, both 
sufficiently simple so as to be realizable in the near fu- 
ture: the linearly inhomogeneous field and the 2D dipole 
field. 



QIDC 



SAW 







^ 




ffc^ 



Moving dot 



>. 




^X 



FIG. 6: Schematic diagram of the spin readout/polarizing de- 
vice based on the Stern-Gerlach effect: SAW propagates from 
left to right transporting a single electron in a moving quan- 
tum dot. QIDC relaxes with a gradient tan{9) to partially 
delocalize the particle during the gate operation; C|(C|) la- 
bels the QIDC receiving electron flux in the spin up(down) 
state of the ay operator; Magnets (of any geometry and polar- 
ity) need to produce a localized, inhomogeneous distribution 
of magnetic field. 



A. Stern-Gerlach gate using a linearly 
inhomogeneous field 

In the first model we will analyze a simple, unidirec- 
tional and linearly inhomogeneous field pointing in the y 
direction: 



B,, 



-f3y, 



(22) 



with Bx = Bz = 0- A wedge-shaped single domain sur- 
face magnet of appropriate dimensions can produce an in- 
homogeneous magnetic field of this form near its center of 
symmetry. A vector potential for this field is A^ = — /3 yz, 
with Ay = Az = 0. This field exerts a spin-dependent 
force in the y direction and it is the simplest field that in- 
duces the Stern-Gerlach effect. Although a z dependence 
enters into the vector potential, by considering motion 
only in the z = plane we may avoid contributions from 
terms involving A in the Hamiltonian, as well as the z 
component of magnetic field. The absence of x in the 
potential immediately allows us to write down harmonic 
oscillator modes for the x dependence. The remaining 
{y, t) dependent part obeys 



1. 



rV^ 



1 



gfiBf^Sy y + VQiDc(y, t)-idt) ip{y, t) = 



(23) 
, and 
Vqidc is given by (|^ . The initial state is again the 
Gaussian ground state of the QIDC. A solution to H23|) 



where Sy is the eigenvalue ±1 of the spinor 



is obtained by a time-dependent Gaussian ansatz 1)11(1 



V'(y, t) = exp [hit) 2/2 + /2(i) y + /3(t)] 



(24) 



In a similar way as before, we solve the system of cou- 
pled ordinary differential equations and derive the time 
dependence of the standard deviation d{t) and the posi- 
tion /i(i) of the probability distribution. 

Let Co = Vq/w'^ and ci = g^iB^Sy. The width of the 
initial Gaussian groundstate wavefunction is determined 
by s = Cq^ . The behavior of the Gaussian solution is 
then characterized by the following time-dependent pa- 
rameters 



d{r) 



f^ir) 



-— ^ [-Aco + a^ cosh(/7 ln(T)/Q) 
-aV^sinh(V7ln(T)/a)] (25) 



y/rci 



2(2a2 + co) 



p/2 _ cosh(V7 hi(T)/a) 



3q 



—- sinh(^ ln(T)/Q)] , 
\/7 



(26) 



where we have introduced ancillary variables a = 
ta,n(6)/w, J — a^ — 4:Co and t = 1 + at. The parameter 
7 is useful in determining when the trigonometric func- 
tions become oscillatory. If the angle 6 is critical such 
that tan(6') = 2y^, then 7 = and the width has an 
especially simple behavior: 



s{tr = - 

a 



l-lnT+-(lnT)^ 



(27) 



This model predicts the width and position of the qubit 
wavefunction after a time t in the linearly inhomogeneous 
field with field gradient /?. In practice we cannot produce 
arbitrary magnetic fields in arbitrary configurations. For 
a weak magnetic field such that cq S> ci 3> ck, 



6{t) - 


" V2V^^ 


= (5(0) W 1 + tan(6l) - 


(28) 


Kr) = 


Cl 

2 Co 


t2- 


~ \Jt cos m r 


■ (29) 



The width increases as ^ \ft and the position moves 
as ~ t^ as expected. We may compare 5(t) with the 
dispersion of a Gaussian in free space 



<5(t) = W<5(0)2 + 



i2 



4 5(0)2 



(30) 



which is linear in t at large times. The broadening of the 
width is suppressed due to the confinement potentials. 
The effect of the Stern-Gerlach can be pictured as fol- 
lows. Consider a qubit in the state | ty) entering the gate. 
The QIDC broadens as the qubit enters the field, allow- 
ing the Stern-Gerlach effect to produce a spin-dependent 
shift in the center of mass towards the channel C| . How- 
ever, the initial localized distribution will deloclize due 



to the QIDC broadening, allowing a small current to be 
detected in the wrong channel C^. Thus this gate is a 
spin polarizing filter with some intrinsic error rate, inde- 
pendent of decoherence effects, shot noise and Johnson 
noise. It remains for us to judge how well it performs, 
and to this end we introduce a fidelity measure of the 
gate. 



B. Fidelity of the Stern-Gerlach gate with linearly 
inhomogeneous field 

As an input to the Stern-Gerlach gate, we consider a 
density matrix of the form 



p(y,2/') = .9(y) 3(2/')* <» cr. 



(31) 



where g(vi) is the Gaussian groundstate of the dot and a 
is the initial spin density matrix. In the inhomogeneous 
field the spatial degrees of freedom couple to the spin 
degrees of freedom, leading to 

P{y,y',t)= Yl 9s{y,t)gs'{y',t)*ass'\s){s'\. (32) 

s,s'=T,i 

After a gate time t, we post-select the state in, say, the 
Ci channel and renormalize: 



dzdz'E^{y,z)gs{z,t)g,,{z',trE^{z\y'), (33) 



The probability yield in the channel is the normalization 
N{t) of the post-selected state. The simplest projection 
kernel projecting into the support of C| is E^{y,y') = 
5{y — y') u{y), where u{y) is the step-function defined in 
(|21|l . We finally define our fidelity as 



dy{]\~p{y.y,t)\]) 



fJtt 

W)Jo 



\g^{y,t)\''dy. 



(34) 

(35) 



The yield Nit) can be calculated without explicit knowl- 
edge of Qi, since the Gaussian solutions are related by 
refiection symmetry: g'\{y,t) = g\_{—y,t). We choose var- 
ious magnetic field gradients ranging from 0.1 T/xm~^ to 
100 T/xm~^ and plot the output fidelity of a maximally- 
mixed input state as a function of the QIDC angle 9 and 
the gate time t (see Fig. [T)). Perfect spin filtering corre- 
sponds to a fidelity F = 1, and the worst-case fidelity is 
F — 1/2. The density plots show how the optimal param- 
eter regions increase with larger magnetic field gradients. 
This model is useful in determining the appropriate gate 
geometry and field gradients required to achieve good fi- 
delities. 

The technology of nanoscale single domain magnets 
is not yet capable of producing large fields in arbitrary 
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FIG. 7: Fidelity F^ in the C| channel for a Stern-Gerlach 
gate operated by a linearly inhomogeneous and unidirectional 
magnetic field in the y direction. White regions => Fj = 1. 9 
is the angle between the QIDC and SAW direction, t is the 
gate time. Input state is a maximally-mixed spin state in the 
groundstate of the SAW quantum dot. Various magnetic field 
gradients have been chosen to show the increase in the optimal 
parameter region (white) in the (f , 9) plane with increasing 
field gradient, imax ~ 100 ns is a crude order of magnitude 
estimate for the GaAs spin relaxation lifetime. 



configurations. The above simple linearly inhomogeneous 
magnetic field is therefore a simple approximation to the 
real magnetic field configuration on any particular device. 
In the next subsection we consider the effect of a different 
(perhaps more realistic) field configuration, namely, the 
field near one of the poles of a dipole magnet. 



C. Stern-Gerlach gate with dipole field 

In this subsection we study the quality of the readout 
obtained from a Stern-Gerlach type gate driven by a sim- 
ple 2D dipole field lying in the 2DEG plane, as shown in 
Fig. |S1 In Appendix ^ we derive the magnetic vector 
potential 



A, 



Px 



{y - df 



(36) 



where /3 is a strength parameter and d is the distance of 
the 2D dipole from the QIDC. 

Figure|51shows that the gradient of the field in the y di- 
rection reverses direction at two points, x = ±d. In order 
to prevent the moving qubit from experiencing opposite 
field gradients, we will restrict the gate time to T = 2 d 




FIG. 8: Vector field plot in the xy (2DEG) plane showing 
the magnetic field due to an infinite string of dipole moments 
running parallel to the z-axis at a distance d from the QIDC. 



so that the qubit moves in the region — T/2 < x < T /2 
in the laboratory frame. The question of how we are to 
shield the extraneous field from the qubit does not have 
a simple solution. It has been suggested that a supercon- 
ducting material such as Niobium could be used to shield 
the magnetic field where it is not neededi^ 

The results of a ID simulation of the wavefunction 
along the x = z — Q direction are shown in Fig. IHl The 
numerical method used was a Crank-Nicholson algorithm 
adapted to the Pauli Equation. The plot shows three 
time-shots of the probability density |V'(2/)^|, normalized 
to peak. The initial state is a Gaussian groundstate with 
spin state | ty). The simulation shows no significant spa- 
tial displacement of the electron density in the y direc- 
tions, in contrast to the behavior exhibited under the lin- 
early inhomogeneous field (c.f. Eq. H26|l ). This is partially 
due to the effect of the Lorentz force generated by the A^ 
term in the Hamiltonian, which contributes an effective 
potential that is highly confining in the y direction. The 
A^ confinement is so strong that the QIDC opening an- 
gle 9 becomes irrelevant. Bohr and Pauli's claim that 
the Lorentz force washes out the Stern-Gerlach effect are 
upheld in this particular scenario. 

Another source of problems for this field configuration 
is that the | ty), | iy) states are not eigenstates of the field 
as they were with the linearly inhomogeneous field. Fig- 
ure ^1 shows the evolution of the Bloch vector in a spin- 
only model of the electron spin qubit. It demonstrates 
how the Bloch vector begins precessing about the x com- 
ponent of the magnetic field, causing the Stern-Gerlach 
force to act in different directions. Although careful tim- 
ing of the gate and possible compensation mechanisms 
could be put in place, this gate is essentially not robust. 
We draw the conclusion that in parameter regimes rele- 
vant to single-electron transport by SAW in GaAs, a field 
of this type will not be suitable for a Stern-Gerlach mea- 
surement gate, and that field unidirectionality is gener- 
ally an important requirement for quantum gates driven 
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FIG. 9: Three time-shots (see legend) showing the spatial 
electron probability density (normalized to peak) . y is the di- 
rection transverse to the QIDC. This is a Stern-Gerlach type 
readout gate using the 2D dipole field (1 Tfj,m~ at x — y = 0) 
shown in Fig. |H| Initial condition is a Gaussian wavepacket 
with spin | ty) • Sideways translation of the probability density 
is suppressed due to both the strongly confining effective po- 
tential arising from the A^ term and the precessional motion 
of the Bloch vector in the non-unidirectional field. 




nuclear spinsi and radio-frequency photons^. 

The models we have presented here give rise to behav- 
ior which was previously neglected in the original pro- 
posal for the quantum computation scheme.';?^ They allow 
us to predict the contribution to decoherence from ex- 
citation into higher orbital states and give their proba- 
bilities of occupation. Any potential decoherence due to 
tunneling between neighboring dots is negligible, and if 
needed the spatial separation of qubits can be increased 
by either introducing higher harmonics of the SAW fun- 
damental frequency or by active gating at the entrance to 
each QIDC. A major source of error is the current fluctu- 
ations in the quantized single-electron SAW current. In 
recent experiments (22iS error in the current quantization 
(including shot and Johnson noise) of less than 0.1% were 
observed. 

The relevant decoherence timescales that underpin our 
proposals is the spin lifetime (Ti and T2) of single elec- 
trons in quantum dots. For n-type bulk semiconduc- 
tors, Ti spin lifetimes in GaAs of 100 ns have been 
reported, which gives a very crude ball-park figured 
An estimate that is more relevant to the SAW electron 
system is the spin relaxation lifetime of a few nanosec- 
onds, which was obtained from spin-resolved micropho- 
toluminescence spectroscopy measurements on photoex- 
cited electronsJ?i However, this is the worst-case scenario 
because the method of electron capture we propose^ pro- 
duces a conduction-band hole that is more extended, 
short-lived and mixed than the photoexcited hole. In- 
deed, recent experiments on static GaAs quantum dots 
report a lower bound on Ti as long as 50 /xs,^^ which 
makes an estimate of 100 ns for the Ti lifetime of SAW 
electrons quite reasonable in the light of current technol- 
ogy. 



FIG. 10: (Color online) Evolution of the Bloch vector in a 2D 
dipole field (1 T/im~^ at x = y = 0): (a^) (long-short dash), 
(ay) (solid), ((Tz) (short dash). Time t — corresponds to the 
starting point x — ~d (see Fig.|HJ. We have used a spin-only 
model of the electron spin qubit to illustrate spin precession 
in the non-unidirectional field. The Bloch vector changes di- 
rection rapidly, which in turn causes the Stern-Gerlach force 
to fluctuate during the operation of the gate. 



by static magnetic fields. 



V. DECOHERENCE EFFECTS 

The major effects of decoherence on this SAW quantum 
computing system have already been considered qualita- 
tively in the original proposal and subsequent papers;^*^ 
and more specifically by a number of authors. These ef- 
fects include: The interaction of the electron qubit with 
other electrons in the 2DEG, surface gates and donor 



VI. SUMMARY 

We have analyzed in detail a proposal^ for imple- 
menting quantum computation on electron spin qubits 
trapped in SAW-QIDC electrostatically-defined dots, us- 
ing only static magnetic gates to perform single-qubit 
initializations, rotations and readouts. Applying the full 
Pauli Hamiltonian, we described the quantum dynamics 
of both the spin and orbital states of the qubit for various 
parameter regimes relevant in SAW single-electron trans- 
port experiments. In the analysis of single-qubit uni- 
tary operations with localized uniform magnetic fields, 
we showed that the effect of the Lorentz force puts upper 
bounds on field strengths and gate times. Moreover, sim- 
ulations showed that field directions normal to the 2DEG 
excite rotational states in the dot. Probabilities of exci- 
tation into higher-energy orbitals were given. In terms 
of feasibility, the models indicate that a field strength of 
about 80 mT will be sufficient to conduct a tt rotation 
in about 1 ns, without compromising the confinement 
properties of the trapped qubit. Since Ti spin lifetimes 
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of microseconds have been reported^^ it is feasible for 
hundreds of single-qubits gates to operate before all co- 
herence is lost in the computation. 

We also studied a device for singie-qubit measurement 
and initialization based on the Stern-Gerlach effect. The 
problematic Lorentz force can be partly suppressed by 
virtue of the geometry of the 2DEG system, and for a 
unidirectional, linearly inhomogeneous field, the correla- 
tion between the spin states and spatial location of the 
qubit leads to a good quantum measurement of its spin. 
For a 2D dipole field, the vector potential has a delete- 
rious effect. Namely, it contributes an effective confin- 
ing potential via the A^ term in the Pauli Hamiltonian 
which suppresses the transverse motion in the y direc- 
tion. Furthermore, the component of magnetic field in 
the X direction causes the Bloch vector to rotate in un- 
desired directions. The latter problem could be over- 
come by using unitary gates prior to the Stern-Gerlach 
gate to correct any spurious rotations, but it is unlikely 
that such a finely-tuned arrangement can be made ro- 
bust. Magnetic fields that are good approximations to a 
unidirectional and linearly inhomogeneous field can pro- 
vide a source of polarized electron qubit states in chan- 
nels Cf and C| with high yield. The very same gate can 
be used to measure the ratio |ap/|/3p for an input qubit 
state \ip) — a\ ly) + f3\ ly). The advantages of this read- 
out method are that the averaging time of the measured 
current can be made long enough to render most noise 
sources (e.g. shot noise, input-referred noise of the cur- 
rent pre-amp) unimportant, and that missing electrons 
from SAW minima do not contribute any error since we 
are measuring the ratio of currents out of the two chan- 
nels. Field gradients in the range 0.1 — 100 Tfim^^ lead 
to gate times that are easily within the range of Ti spin 
lifetimes in GaAs. These parameter regimes will need to 
be probed by experiment if we are to make progress to- 
wards feasible quantum computation with nanomagnets 
in the SAW quantum computer. 
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APPENDIX A: DERIVATION OF THE 2D 
DIPOLE FIELD VECTOR POTENTIAL 



We adopt cylindrical coordinates (r, 0, y) about the y 
axis, with r — \/ x^ + z^ . We consider the field due to 
an infinite string of dipoles passing parallel to the z-axis 
through y = d as shown in Fig. |H1 The resulting mag- 
netic field has planar symmetry with respect to the xy 
(2DEG) plane. The QIDC runs along the x axis. A 
vector potential for this field can be constructed as fol- 
lows. A magnetic dipole at y = z = pointing in the +y 
direction has a vector potential — 



A'4r,y) 



jiQiriBdz 



An 



(r2 +y2)3/2' 



(Al) 



where A'^. — A' — and uib is the dipole moment per 
unit length in the z direction. We now displace this po- 
tential by d in the y direction and integrate over all z, 
obtaining 



A, = 



/3x 



{y-df 



(A2) 



with At 



A,„ 



and strength parameter /3 



fiom'g/2Tr. The magnetic field can be derived by applying 
the curl in cylindrical coordinates: Br — —r~^dy{rAci,) 
^■^dr{rA^). Figure IHl shows the vector field 



and B 



y 



plot of the resulting magnetic field, which falls off as 
for large r. 
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